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Abstract 



Using the Landauer formulation of transport theory and tight binding models 
of the electronic structure, we study electron transport through atomic wires 
that form ID constrictions between pairs of metallic nano-contacts. Our re- 
sults are interpreted in terms of electron standing waves formed in the atomic 
wires due to interference of electron waves reflected at the ends of the atomic 
constrictions. We explore the influence of the chemistry of the atomic wire- 
metal contact interfaces on these standing waves and the associated transport 
resonances by considering two types of atomic wires: gold wires attached to 
gold contacts and carbon wires attached to gold contacts. We find that the 
conductance of the gold wires is roughly IGq = 2e 2 /h for the wire lengths 
studied, in agreement with experiments. By contrast, for the carbon wires 
the conductance is found to oscillate strongly as the number of atoms in the 
wire varies, the odd numbered chains being more conductive than the even 
numbered ones, in agreement with previous theoretical work that was based 
on a different model of the carbon wire and metal contacts. 
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I. INTRODUCTION 



The ability to manipulate matter on the atomic scale has improved greatly over the past 
decade. The STM has facilitated much of this progress with its ability to positionally control 
atoms on surfaces.0~i Experiments are now able now to move beyond the manipulation stage 
and measure the electronic properties of these atomic scale systems. One system of particular 
interest is an atomic chain. When an atomic chain is used to bridge two contacts an atomic 
wire is formed. Recent STM experiments have managed to measure the tunneling current 
through atomic wires, which have been made by positioning several atoms to form a chain 
between the STM tip and the substrate. i@ More recently, mechanically controlled break 
junctions have been used to form a gold chain between two gold contacts. The conductance 
of these gold atomic wires was measured for different chain lengths J It may also be possible 
to form carbon atomic wires by stretching a single walled nanotube.i Current work is focusing 
on connecting carbon wires between two metalized AFM tips.0 Early theoretical work on 
atomic wires studied the electron transmission through molecular chains with a single orbital 
per site with nearest neighbour hopping.0 Recent theoretical studies of these atomic wire 
systems based on a density functional approach have also revealed interesting phenomena, 
one of them being the oscillation of conductance versus the length of a carbon wire. El 

The atomic wire is interesting from a transport point of view because it confines electrons 
to propagate only in one dimension. It is also on a length scale that is of the order of 
the wavelength of the electrons which are incident from the contacts. This means that 
the quantum mechanical wave nature of the electrons will be important. When waves are 
confined to propagate in a ID constriction standing waves are formed. These arise because 
of back reflections at the boundaries between the contacts and the constriction. Electrons 
propagating in semiconductor waveguides have been predicted to form standing waves along 
the conduction channel Jllll Experimentally observing standing wave formation in these 
systems has been difficult because these structures are formed electro-statically, which results 
in ID channels that are smoothly tapered at their ends minimizing reflections there. BiWe 
show that electron standing waves should form in atomic wires and investigate the role of the 
chemistry of the wire-contact interface (which has no analog in the semiconductor devices) 
on their behavior. Although these standing waves have not been studied experimentally at 
the present time, atomic wire systems may prove to be better candidates for their observation 
since the interfaces between a ID wire and the contacts cannot be smooth on the atomic 
scale so that reflections at these interfaces should be significant. 

We study electron transport through an atomic wire by calculating the wavefunction for 
an electron propagating along the wire. This is done using a tight binding formalism that 
treats both the chain and metallic leads atomisticallyJil The present work is complemen- 
tary to other theoretical treatments that have been based on tight binding models of the 
atomic wirefflJll or have considered jellium contacts with the wire modelled using pseudo- 
potentials.il! We apply our technique to two different atomic wire systems. The first system 
is a homogeneous atomic wire which is made using a gold chain that bridges two gold con- 
tacts. We find that our results are consistent with the interpretation that standing waves 
are formed along the gold wire. The standing wave patterns that are formed are confined 
within the ID channel that has the length of the gold chain. The number of standing wave 
resonances that occurs increases with the length of the chain. However, the conductance 
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of gold chains consisting of different numbers of atoms is found to vary only slightly as the 
number of atoms changes and remains near unity as has been reported experiment allyi 

The second system we consider is a heterogeneous atomic wire. It is a carbon atom 
chain attached to two gold contacts. Again our results are consistent with a description 
involving the formation of standing waves along the carbon chain. However, in this case the 
standing wave patterns formed extend significantly beyond the ends of the chain into the 
gold contacts. Thus in this heterogeneous system the chain is effectively lengthened due to 
the chemistry between the chain and the contacts. This gives rise to a greater number of 
standing wave resonances. The conductance is also found to oscillate with the length of the 
chain which is in agreement with what was found in other work .111 The standing wave model 
that we describe provides an intuitive physical explanation of this predicted oscillation. 

In Sec. I we present the theoretical approach used to calculate the conductance together 
with a description of standing wave formation in an atomic wire. We then present our 
calculations for gold atomic wires, a homogeneous system, in Sec. II. In Sec. Ill we present 
calculations for carbon atomic wires. 



II. ELECTRON TRANSPORT IN ATOMIC WIRES: NORMAL MODES 

An atomic wire is composed of a chain of atoms that forms a ID constriction between two 
metallic contacts which act as the source and the drain. Electrons incident as Bloch waves 
from the source contact scatter through the ID constriction and enter the drain. Landauer 
theory relates the conductance G of the atomic wire to the multi-channel transmission 
probability T for an electron at the Fermi level of the source to scatter from source to drain 
via G = ^T.§ 

h 

We calculate the multi-channel transmission probability T by solving for the t matrix of 
the atomic wire system. For a given energy E, an electron can propagate in one of many 
electronic modes in both the source and the drain. The t matrix connects the modes (or 
channels) in the source to the channels in the drain for a given energy. We use a tight- 
binding treatment of Schrodinger's equation to set up a system of linear equations to solve 
for t. The source and drain are treated as ID leads with multi-atom unit cells. The wire is 
attached to atomic tips which are then attached to the leads. An schematic diagram of an 
atomic wire system is shown in Fig. |V|. 

The first step in the transport calculation is to impose boundary conditions on the 
wavefunction \^ a ) which describes an electron incident in the a th channel with energy E. 
In the source (or left lead (L)) consists of the incident rightward propagating (+) a th 
Bloch mode along with a sum over reflected leftward propagating (-) and evanescent 

modes On the atomic wire (A), the wavefunction is a sum over the atomic orbitals 

\j) that exist on the wire. In the drain (or right lead (R)) the wavefunction is a sum over the 
transmitted rightward propagating and evanescent modes |<&" )r. Thus the wavefunction is 
given by, 

|* a ) = |* a )L+|* >A+|* a >fl (1) 

where 
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\y a ) L = i$:)L+E r «',«i $ -)i (2) 

a' 

\* Q )A = Y. d j\j) (3) 

3 

l* a >K = E*«*',«K>* ( 4 ) 
a' 

Using Schrodinger's equation if |\I/ a ) = E\ty a ), where H is the extended Hiickel tight binding 
HamiltoniarS of the system, yields a system of equations which we solve to find the r's, ef s 
and fs. With the t matrix thus determined, the transmission probability T is then found 
using 

n^) = EE v i \HA 2 (5) 

where the sum over a is over the rightward propagating modes in the left lead and the sum 
over (3 is over rightward propagating modes in the right lead. The velocity ratio appears 
since the velocities of modes in the left and right leads may be different. The conductance 
is then found by using the Landauer formula given above for the transmission probability 
evaluated at the Fermi energy Ef of the source. This method and the procedure used for 
handling the non-orthogonality between the different tight-binding orbitals is described in 
more detail elsewhere .£30 

The above tight-binding model is completely general and is applicable to many types 
of systems, not just atomic wires. We now proceed to examine the physical behaviour the 
above model should exhibit when applied to a ID constriction such as an atomic wire. The 
physical interpretation that we will use to describe the results of our transport calculations 
that follow is in terms of electronic standing waves formed within the constriction. 

When waves (be they light, sound or electron) are confined to propagate within a finite 
one dimensional channel, standing waves are formed. They are formed due to the interference 
of waves back reflected from the ends of the pipe where there is a large impedance mismatch. 
These standing waves (also known as normal modes) have an associated wave vector given 
by k = nir/L where n — 1, 2, 3, . . . labels the normal mode and L is the effective length of the 
channel. Within the tight binding framework, for a ID channel composed of a finite chain of 
discrete atoms equally spaced by a distance a, the number of normal modes becomes finite. 
The reduced wavevector of the wave is y = ka = rnra/L and has the range < y < tt. Thus 
na must be less than or equal to L which restricts the values that n can take. Thus a finite 
chain supports only a finite number of standing waves. 

In regards to electron transport, when an electron is incident from the source with an 
energy that corresponds to the energy of one of the normal modes, it will form a standing 
wave along the chain. This will correspond to maximal transmission through the atomic wire. 
Because the chain supports a finite number of normal modes the transmission probability 
will have a discrete number of resonances, each one corresponding to one of the normal 
modes. 

The effective length L of the atomic wire determines the number of resonances in the 
transmission probability. Intuitively this length should approximate the contact-to-contact 
distance. However, an effective length which may be shorter or longer than the contact-to- 
contact spacing will occur in practice. This depends on the chemical and geometrical nature 
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of the interface between the atomic chain and the metallic contacts. The chemistry may serve 
to effectively sharpen or broaden the constriction which will affect the effective length. For 
a homogeneous atomic wire system where the contacts and the chain are made of the same 
atoms, the length L of the wire should be approximately equal to the contact-to-contact 
distance. In this situation the confining potential is that created by the geometry of the 
constriction and hence the backscattering should occur at the interface between the contacts 
and the chain where the impedance mismatch is high. For heterogeneous systems, the length 
L may be different from the contact-to-contact length. This is because the chemistry between 
the chain and the contact is significantly different than that within the bulk of the contact. 
Thus the confining potential due to the geometry is now altered. An interface region is 
created which modifies the confining potential and this changes the effective length L of the 
chain. So in heterogeneous wires one may expect there to be more or fewer resonances than 
in the homogeneous case. 

III. HOMOGENEOUS WIRE: GOLD WIRE WITH GOLD CONTACTS 

The first system we consider is a homogeneous atomic wire of gold. An atomistic diagram 
of one of the wires is shown in Fig. |V[ The atomic chain is made up of gold atoms with a 
spacing of 2.4 A. The chain is bonded over two gold tips, each consisting of 13 gold atoms 
and both oriented in the (100) direction. The perpendicular bonding distance between the 
chain and a tip is taken to be 2.0 A. Each tip is then bonded to a layer of 16 gold atoms 
which forms the unit cell for the contact. The gold atoms along the chain and the four gold 
atoms in each tip to which the chain is attached are modeled using the 6 s6 p5 d orbitals of 
gold. The remaining gold atoms in the tips and those in the leads are modeled using a 6s 
orbital. We consider chains ranging in length from 3 to 7 gold atoms. 

The transmission probabilities for the five different chain lengths are shown in Fig. |V]. 
As expected there is a different number of resonances for each chain length. For the energies 
shown, only a 6s mode is supported by the chain and hence the maximum value that the 
transmission probability can take is 1 (ignoring spin). If there were more modes which could 
propagate along the chain the multi-channel Landauer transmission probability could be 
greater than 1. Since the gold atomic wires are a homogeneous system the length corresponds 
approximately to the contact-to-contact spacing which is roughly (N + l)a where N is the 
number of atoms forming the chain and a = 2.4A. So an N atom chain should display N 
resonances and this is seen. (One might argue that there should be N + 1 resonances, but 
the wavefunction that corresponds to this n value has a node on each atomic site and hence 
is unphysical). The N or last resonance of each chain (which is not shown) has an energy 
that lies between -6.5 eV and -5.5 eV and this overlaps some 6 p modes which also contribute 
to the transmission and obscure this resonance. 

At each resonance the electron wavefunction is a standing wave along the 6 s orbital 
backbone of the chain. The squared modulus should behave as sin 2 (nll(x — xl)/L) where 
x is the position coordinate along the wire (x = marks the first chain atom), xl is the 
leftmost position where the wavefunction is zero (this may extend beyond the chain into 
the tips) and n labels the mode. The wavefunction (calculated using our tight-binding 
model) for the 5 atom chain is plotted for the n = 1 and n = 4 resonance together with the 
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corresponding ideal sinusoidal standing waveform in Fig. 0a. The graphs clearly indicate 
the formation of standing waves along the chain. The length L which best fits the data is 
around 13A which is just slightly shorter than the 13.6Acontact-to-contact spacing for this 
chain. In Fig. [V|a the behaviour of the wavefunction in the gold tips for the n = 3 mode is 
depicted for the 5 atom chain. The atoms in the tips fall outside the length of the "pipe" 
and thus are not part of the standing wave formation within the pipe. The wavefunction 
is consistently larger within the left tip than on the right as is to be expected for a wave 
incident on the wire from the left contact. (Note that the squared wavefunction that is 
depicted is a qualitative result and should not be compared to the transmission results. 
What is shown is the result of adding together the squared amplitude of the wavefunction 
for each incident channel. Thus it is the "average" wavefunction along the atomic chain.) 

The conductance of these chains is given by the Landauer formula, G = GqT. The 
transmission probability that enters the Landauer expression is evaluated at the source 
Fermi energy. For our gold leads the Fermi energy is approximately -10 eV. Fig. |V|a shows 
the conductance as a function of the number of atoms in the wire at this Fermi energy. In 
qualitative agreement with experimentll the conductance is nearly unity and varies between 
0.81 G and 0.99 G a . The fluctuations are relatively small because the peak to valley ratio 
of the transmission probability near the Fermi energy is near unity for the different gold 
chains. 

IV. HETEROGENEOUS WIRE: CARBON WIRE RESULTS 

The second structure that we consider is a heterogeneous system consisting of a carbon 
atomic wire and gold contacts. The atomic system is similar to that depicted in Fig. |V]. 
The carbon chain spacing is taken to be 1.32A. Again it is bonded over 13-atom gold tips 
oriented in the (100) direction. The perpendicular bond distance between the chain and a 
tip was chosen to be lA. The four gold atoms that bond to the carbon wire on each tip were 
simulated using their 6 s6p5 d orbitals while the other gold atoms were simulated using just 
their 6 s orbital. (We remark that the results obtained modelling all of the gold atoms in this 
system with only 6s orbitals were qualitatively similar.) The carbon atoms were simulated 
using their 2 s and 2p orbitals. Chains consisting of 4 to 10 carbon atoms were simulated. 

The transmission probabilities for chains numbering 4 through 8 are shown in Fig. [V]. 
The magnitude is greater than that found in the gold chains with a maximum value of near 
2.0. This is because for these energies the electrons are propagating along the tt backbone 
of the carbon wire. If the wire is oriented along the z direction there are two independent 
modes (p x and p y ) along which the electrons can propagate. Thus the transmission can be 
a maximum of 2.0 (ignoring spin). In those regions where it exceeds 2.0, an extra a mode 
is present. The number of resonances for a given numbered carbon chain is also greater 
than that of the corresponding gold chain. For heterogeneous systems the chemistry at the 
interface between the atomic wire and the metallic contacts effectively alters the length of 
the ID conducting channel. For example the 4 atom carbon chain displays 6 resonances 
which shows that the interaction is such that it is effectively lengthened and behaves like a 
6 atom chain. 

To determine what the effective length for one of these carbon chains is we examine the 
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electron's wavefunction. The wavefunctions for two different resonances of the 5 atom carbon 
chain are plotted in Fig. ^b. The corresponding normal modes are also plotted where L 
was chosen so as to give a best fit to all the curves. The value of L was determined to 
be around lO.OAwhich is significantly longer than the contact-to-contact distance of 7.28A. 
From this we can estimate that the interface region extends about 1.36Ainto the gold tip. 
So in this system the standing waves are forming from within the gold tips rather than at 
their surfaces as was the case for the gold atom wires. The behaviour of the wavefunction 
in the tips is depicted in Fig. |V|b. It is seen that the standing wave continues onto the 
first layer of each tip. However, the wavefunction on the second atomic layer of each tip lies 
outside the length L and resembles the behavior of the wavefunction within the tips of the 
gold atomic wire. It is seen that the wavefunction on the left tip is greater than that on the 
right. Again, we emphasize that we are using the "average" wavefunction to perform this 
analysis. 

The conductance of these carbon chains was found using the same source Fermi energy 
of -10 eV. Conductance oscillations are found for this choice as the number of carbon atoms 
in the chain is varied, with odd numbered chains possessing a higher conductance than even 
ones (shown in Fig. 0b). This is consistent with other theoretical work which has also found 
conductance oscillations as a function of carbon wire length although the carbon wire and 
contacts were modelled in a different way in that work.0 The oscillation has quite a strong 
signature with the conductance varying between 0.9 G Q and 1.8 G Q . For the odd numbered 
chains the conductance is large because the Fermi energy corresponds to a standing wave 
resonance, whereas the even numbered chains are off resonance. The standing wave inter- 
pretation seems to be a reasonable physical interpretation for the oscillating conductance 
found here and in a previous report .E3 



V. CONCLUSIONS 

We have shown that electron transport through an atomic wire is consistent with the 
interpretation that electron standing waves are formed in the atomic wire which bridges 
metallic contacts. They arise because of the interference of reflected electron waves at the 
interface between the atomic chain and the metallic tips. For the gold tips the effective 
length of the constriction was roughly the contact-to-contact distance. The conductance of 
the gold wires was also found to be near unity for atomic wires ranging from 3 to 7 gold atoms, 
consistent with experiment. The behavior in the carbon wires was different. The effective 
length of the constriction was significantly larger signalling that the formation of standing 
waves is sensitive to the chemistry of the atomic wire-contact interface. The conductance 
was also found to oscillate as a function of chain length and this can be understood in terms 
of the different standing wave modes that are created in chains of differing length. The 
observation of electron standing waves may be facilitated by using atomic wires (as opposed 
to semiconductor systems) since the constriction can not be smooth, which will enhance the 
reflections. 

This work was supported by NSERC. 
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FIG. 3. (a) Average squared wavefunctions (diamonds) for the n = 1, 4 transmission probability 
resonances of the 5 atom gold chain. The corresponding normal modes, sin 2 (n7r(x — xl)/L), with 
n = 1,4 are also plotted. The variable X measures the length along the chain and is in A. (b) 
Average squared wavefunctions for the n = 3, 4 transmission probability resonances of the 5 atom 
carbon chain. The corresponding normal modes are shown for n = 3,4. 
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FIG. 6. (a) Conductance of N gold atom chains, where N = 3, . . . , 7. (1 G = 2e 2 /h). The 
Fermi energy of the contacts was chosen to be -10.0 eV. (b) Conductance of N carbon atom chains 
bonded to gold leads, where N = 4, . . . , 10. The Fermi energy of the contacts was chosen to be 
-10.0 eV. 
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